Import Packages
library(ggplot2)
library(plotly)
library(data.table)
Import Data and look at first couple records
setwd("F:/kaggle/Mercedes-Benz")
sample_submission <- fread("data/sample_submission.csv")
train <- fread("data/train.csv")
test <- fread("data/test.csv")
train[1:2,]
# 1 - ID
# 2 - Response Varuable (Time in Seconds on Test Stand)
# 3-10 String Option Codes
# 11-385 0/1 values based on option codes
# X7 and X9 are missinf for some reason.
Calculate Mean and SD for Response Variable (Y). This is the time (sec) for a vehicle on the MB test station.
mean(train[,y])
[1] 100.6693
sd(train[,y])
[1] 12.67938
Create plots for the Response Variable. Sorted by ID.
plot_ly(y=train$y, type="scatter")
Create plots for the Response Variable. Sorted by the Response Variable.
plot_ly(y=train[order(y),]$y, type="scatter")
Create histogram for the Response Variable.
plot_ly(x=train$y, type="histogram")
Create histogram for the LOG(Response Variable).
plot_ly(x=log(train$y), type="histogram")
Frequency tables for the first 10 variables (All with String Codes)
cbind(train[,.N,by=X0],Mean_Y=round(train[,mean(y),by=X0]$V1,1))[order(-N)]
cbind(train[,.N,by=X1],Mean_Y=round(train[,mean(y),by=X1]$V1,1))[order(-N)]
cbind(train[,.N,by=X2],Mean_Y=round(train[,mean(y),by=X2]$V1,1))[order(-N)]
cbind(train[,.N,by=X3],Mean_Y=round(train[,mean(y),by=X3]$V1,1))[order(-N)]
cbind(train[,.N,by=X4],Mean_Y=round(train[,mean(y),by=X4]$V1,1))[order(-N)]
cbind(train[,.N,by=X5],Mean_Y=round(train[,mean(y),by=X5]$V1,1))[order(-N)]
cbind(train[,.N,by=X6],Mean_Y=round(train[,mean(y),by=X6]$V1,1))[order(-N)]
cbind(train[,.N,by=X8],Mean_Y=round(train[,mean(y),by=X8]$V1,1))[order(-N)]#where is 7
#where is 9
Frequency plots for the first 10 variables (All with String Codes)
# column_names_for_option_plots_string_codes <- colnames(train)[3:5] #3:12
# for(i in column_names_for_option_plots_string_codes){
# x_values <- as.data.frame(train[,.N,by=i][order(-N)])[,i]
# y_values <- as.data.frame(train[,.N,by=i][order(-N)])[,'N']
# plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar")
##Couldn't get loop plotting to work => https://github.com/ropensci/plotly/issues/273
#Plotting manually
column_names_for_option_plots_string_codes <- colnames(train)[3:10] #3:10
i<-1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
i<-i+1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
i<-i+1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
i<-i+1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
i<-i+1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
i<-i+1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
i<-i+1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
i<-i+1
current_colname <- column_names_for_option_plots_string_codes[i]
x_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,current_colname]
y_values <- as.data.frame(train[,.N,by=current_colname][order(-N)])[,'N']
plot_ly(x=as.list(x_values),y=as.list(y_values), type="bar") %>% layout(title=current_colname)
Frequency table for the remaining variables (All with 0/1 Coding)
column_names_for_summary <- colnames(train)[11:length(colnames(train))]
Variable <- c()
N_0 <- c()
N_1 <- c()
Mean_0 <- c()
Mean_1 <- c()
for(i in column_names_for_summary){
j<-1
Variable <- c(Variable,i)
N_0 <- c(N_0,train[,.N,by=i][1]$N)
N_1 <- c(N_1,train[,.N,by=i][2]$N)
Mean_0 <- c(Mean_0,round(train[,mean(y),by=i][1]$V1,1))
Mean_1 <- c(Mean_1,round(train[,mean(y),by=i][2]$V1,1))
j <- j+1
}
# N_0[is.na(N_0)==TRUE] <- 0
# N_1[is.na(N_1)==TRUE] <- 0
# Mean_0[is.na(Mean_0)==TRUE] <- 0
# Mean_1[is.na(Mean_1)==TRUE] <- 0
summary_results <- as.data.frame(cbind(Variable,N_0=N_0,N_1,Mean_0,Mean_1), stringsAsFactors = FALSE)
summary_results$N_0 <- as.integer(summary_results$N_0)
summary_results$N_1 <- as.integer(summary_results$N_1)
summary_results$Mean_0 <- as.numeric(summary_results$Mean_0)
summary_results$Mean_1 <- as.numeric(summary_results$Mean_1)
summary_results$Delta_Mean <- summary_results$Mean_1 - summary_results$Mean_0
summary_results
After Looking at Kaggle, checked for Duplicate fileds - Added Data to Frequency table for the remaining variables (All with 0/1 Coding).
train_2 <- train[, !duplicated(t(train))] #remove duplicated fields ... from raddar@Kaggle => https://www.kaggle.com/c/mercedes-benz-greener-manufacturing/discussion/34006
duplicate_column <- c()
for(i in duplicated(t(train))){
duplicate_column <- c(duplicate_column,i)
}
duplicate_column <- duplicate_column[11:length(duplicate_column)]
summary_results$duplicate_column <- duplicate_column
summary_results_decreasing <- summary_results[order(-N_1),]
summary_results
summary_results_decreasing
List of Duplicated Columns
summary_results_decreasing[(summary_results_decreasing$duplicate_column==TRUE),]
Create Response Plots for all Binary Options
train_df <- as.data.frame(train)
train_df_dedupe <- train_df[, !duplicated(t(train_df))] #remove duplicated fields ... from raddar@Kaggle => https://www.kaggle.com/c/mercedes-benz-greener-manufacturing/discussion/34006
column_names_for_binary_plots <- colnames(train_df_dedupe)[11:length(colnames(train_df_dedupe))]
l <- htmltools::tagList()
for (i in 1:length(column_names_for_binary_plots)) { #
plot_variable <- column_names_for_binary_plots[i]
l[[i]] <- plot_ly(x=train_df_dedupe[,plot_variable],y=train_df_dedupe$y, type="box", height=300) %>% layout(title=plot_variable)
}
l
LS0tDQp0aXRsZTogIk1lcmNlZGVzLUJlbnogRURBIE5vdGVib29rIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQphdXRob3I6IEplZmYgSGVkYmVyZw0KZGF0ZTogNS1KdW5lLTIwMTcgICAgDQotLS0NCjxicj4NCg0KIyMjIyBJbXBvcnQgUGFja2FnZXMNCmBgYHtyLG1lc3NhZ2U9RkFMU0Usd2FybmluZz1GQUxTRX0NCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkocGxvdGx5KQ0KbGlicmFyeShkYXRhLnRhYmxlKQ0KYGBgDQoNCjxicj4NCg0KIyMjIyAgSW1wb3J0IERhdGEgYW5kIGxvb2sgYXQgZmlyc3QgY291cGxlIHJlY29yZHMNCmBgYHtyfQ0Kc2V0d2QoIkY6L2thZ2dsZS9NZXJjZWRlcy1CZW56IikNCg0Kc2FtcGxlX3N1Ym1pc3Npb24gPC0gZnJlYWQoImRhdGEvc2FtcGxlX3N1Ym1pc3Npb24uY3N2IikNCnRyYWluIDwtIGZyZWFkKCJkYXRhL3RyYWluLmNzdiIpDQp0ZXN0IDwtIGZyZWFkKCJkYXRhL3Rlc3QuY3N2IikNCg0KdHJhaW5bMToyLF0NCg0KIyAxIC0gSUQNCiMgMiAtIFJlc3BvbnNlIFZhcnVhYmxlIChUaW1lIGluIFNlY29uZHMgb24gVGVzdCBTdGFuZCkNCiMgMy0xMCAgU3RyaW5nIE9wdGlvbiBDb2Rlcw0KIyAxMS0zODUgMC8xIHZhbHVlcyBiYXNlZCBvbiBvcHRpb24gY29kZXMNCg0KIyBYNyBhbmQgWDkgYXJlIG1pc3NpbmYgZm9yIHNvbWUgcmVhc29uLg0KYGBgDQo8YnI+DQoNCiMjIyMgIENhbGN1bGF0ZSBNZWFuIGFuZCBTRCBmb3IgUmVzcG9uc2UgVmFyaWFibGUgKFkpLiAgVGhpcyBpcyB0aGUgdGltZSAoc2VjKSBmb3IgYSB2ZWhpY2xlIG9uIHRoZSBNQiB0ZXN0IHN0YXRpb24uDQpgYGB7cn0NCm1lYW4odHJhaW5bLHldKQ0Kc2QodHJhaW5bLHldKQ0KYGBgDQo8YnI+DQoNCiMjIyMgIENyZWF0ZSBwbG90cyBmb3IgdGhlIFJlc3BvbnNlIFZhcmlhYmxlLiAgU29ydGVkIGJ5IElELg0KYGBge3IsIG1lc3NhZ2U9RkFMU0V9DQpwbG90X2x5KHk9dHJhaW4keSwgdHlwZT0ic2NhdHRlciIpDQpgYGANCjxicj4NCg0KIyMjIyAgQ3JlYXRlIHBsb3RzIGZvciB0aGUgUmVzcG9uc2UgVmFyaWFibGUuICBTb3J0ZWQgYnkgdGhlIFJlc3BvbnNlIFZhcmlhYmxlLg0KYGBge3IsIG1lc3NhZ2U9RkFMU0V9DQpwbG90X2x5KHk9dHJhaW5bb3JkZXIoeSksXSR5LCB0eXBlPSJzY2F0dGVyIikNCmBgYA0KPGJyPg0KDQojIyMjICBDcmVhdGUgaGlzdG9ncmFtIGZvciB0aGUgUmVzcG9uc2UgVmFyaWFibGUuDQpgYGB7cn0NCiBwbG90X2x5KHg9dHJhaW4keSwgdHlwZT0iaGlzdG9ncmFtIikNCmBgYA0KPGJyPg0KDQojIyMjICBDcmVhdGUgaGlzdG9ncmFtIGZvciB0aGUgTE9HKFJlc3BvbnNlIFZhcmlhYmxlKS4NCmBgYHtyfQ0KIHBsb3RfbHkoeD1sb2codHJhaW4keSksIHR5cGU9Imhpc3RvZ3JhbSIpDQpgYGANCjxicj4NCg0KIyMjIyAgRnJlcXVlbmN5IHRhYmxlcyBmb3IgdGhlIGZpcnN0IDEwIHZhcmlhYmxlcyAoQWxsIHdpdGggU3RyaW5nIENvZGVzKQ0KYGBge3J9DQpjYmluZCh0cmFpblssLk4sYnk9WDBdLE1lYW5fWT1yb3VuZCh0cmFpblssbWVhbih5KSxieT1YMF0kVjEsMSkpW29yZGVyKC1OKV0NCmNiaW5kKHRyYWluWywuTixieT1YMV0sTWVhbl9ZPXJvdW5kKHRyYWluWyxtZWFuKHkpLGJ5PVgxXSRWMSwxKSlbb3JkZXIoLU4pXQ0KY2JpbmQodHJhaW5bLC5OLGJ5PVgyXSxNZWFuX1k9cm91bmQodHJhaW5bLG1lYW4oeSksYnk9WDJdJFYxLDEpKVtvcmRlcigtTildDQpjYmluZCh0cmFpblssLk4sYnk9WDNdLE1lYW5fWT1yb3VuZCh0cmFpblssbWVhbih5KSxieT1YM10kVjEsMSkpW29yZGVyKC1OKV0NCmNiaW5kKHRyYWluWywuTixieT1YNF0sTWVhbl9ZPXJvdW5kKHRyYWluWyxtZWFuKHkpLGJ5PVg0XSRWMSwxKSlbb3JkZXIoLU4pXQ0KY2JpbmQodHJhaW5bLC5OLGJ5PVg1XSxNZWFuX1k9cm91bmQodHJhaW5bLG1lYW4oeSksYnk9WDVdJFYxLDEpKVtvcmRlcigtTildDQpjYmluZCh0cmFpblssLk4sYnk9WDZdLE1lYW5fWT1yb3VuZCh0cmFpblssbWVhbih5KSxieT1YNl0kVjEsMSkpW29yZGVyKC1OKV0NCmNiaW5kKHRyYWluWywuTixieT1YOF0sTWVhbl9ZPXJvdW5kKHRyYWluWyxtZWFuKHkpLGJ5PVg4XSRWMSwxKSlbb3JkZXIoLU4pXSN3aGVyZSBpcyA3DQojd2hlcmUgaXMgOQ0KDQpgYGANCjxicj4NCg0KIyMjIyAgRnJlcXVlbmN5IHBsb3RzIGZvciB0aGUgZmlyc3QgMTAgdmFyaWFibGVzIChBbGwgd2l0aCBTdHJpbmcgQ29kZXMpDQpgYGB7cn0NCiMgY29sdW1uX25hbWVzX2Zvcl9vcHRpb25fcGxvdHNfc3RyaW5nX2NvZGVzIDwtIGNvbG5hbWVzKHRyYWluKVszOjVdICMzOjEyDQojIGZvcihpIGluIGNvbHVtbl9uYW1lc19mb3Jfb3B0aW9uX3Bsb3RzX3N0cmluZ19jb2Rlcyl7DQojICAgeF92YWx1ZXMgPC0gYXMuZGF0YS5mcmFtZSh0cmFpblssLk4sYnk9aV1bb3JkZXIoLU4pXSlbLGldDQojICAgeV92YWx1ZXMgPC0gYXMuZGF0YS5mcmFtZSh0cmFpblssLk4sYnk9aV1bb3JkZXIoLU4pXSlbLCdOJ10NCiMgICBwbG90X2x5KHg9YXMubGlzdCh4X3ZhbHVlcykseT1hcy5saXN0KHlfdmFsdWVzKSwgdHlwZT0iYmFyIikNCiMjQ291bGRuJ3QgZ2V0IGxvb3AgcGxvdHRpbmcgdG8gd29yayA9PiBodHRwczovL2dpdGh1Yi5jb20vcm9wZW5zY2kvcGxvdGx5L2lzc3Vlcy8yNzMNCg0KI1Bsb3R0aW5nIG1hbnVhbGx5DQpjb2x1bW5fbmFtZXNfZm9yX29wdGlvbl9wbG90c19zdHJpbmdfY29kZXMgPC0gY29sbmFtZXModHJhaW4pWzM6MTBdICMzOjEwDQppPC0xDQpjdXJyZW50X2NvbG5hbWUgPC0gY29sdW1uX25hbWVzX2Zvcl9vcHRpb25fcGxvdHNfc3RyaW5nX2NvZGVzW2ldDQp4X3ZhbHVlcyA8LSBhcy5kYXRhLmZyYW1lKHRyYWluWywuTixieT1jdXJyZW50X2NvbG5hbWVdW29yZGVyKC1OKV0pWyxjdXJyZW50X2NvbG5hbWVdDQp5X3ZhbHVlcyA8LSBhcy5kYXRhLmZyYW1lKHRyYWluWywuTixieT1jdXJyZW50X2NvbG5hbWVdW29yZGVyKC1OKV0pWywnTiddDQpwbG90X2x5KHg9YXMubGlzdCh4X3ZhbHVlcykseT1hcy5saXN0KHlfdmFsdWVzKSwgdHlwZT0iYmFyIikgJT4lIGxheW91dCh0aXRsZT1jdXJyZW50X2NvbG5hbWUpDQppPC1pKzENCmN1cnJlbnRfY29sbmFtZSA8LSBjb2x1bW5fbmFtZXNfZm9yX29wdGlvbl9wbG90c19zdHJpbmdfY29kZXNbaV0NCnhfdmFsdWVzIDwtIGFzLmRhdGEuZnJhbWUodHJhaW5bLC5OLGJ5PWN1cnJlbnRfY29sbmFtZV1bb3JkZXIoLU4pXSlbLGN1cnJlbnRfY29sbmFtZV0NCnlfdmFsdWVzIDwtIGFzLmRhdGEuZnJhbWUodHJhaW5bLC5OLGJ5PWN1cnJlbnRfY29sbmFtZV1bb3JkZXIoLU4pXSlbLCdOJ10NCnBsb3RfbHkoeD1hcy5saXN0KHhfdmFsdWVzKSx5PWFzLmxpc3QoeV92YWx1ZXMpLCB0eXBlPSJiYXIiKSAlPiUgbGF5b3V0KHRpdGxlPWN1cnJlbnRfY29sbmFtZSkNCmk8LWkrMQ0KY3VycmVudF9jb2xuYW1lIDwtIGNvbHVtbl9uYW1lc19mb3Jfb3B0aW9uX3Bsb3RzX3N0cmluZ19jb2Rlc1tpXQ0KeF92YWx1ZXMgPC0gYXMuZGF0YS5mcmFtZSh0cmFpblssLk4sYnk9Y3VycmVudF9jb2xuYW1lXVtvcmRlcigtTildKVssY3VycmVudF9jb2xuYW1lXQ0KeV92YWx1ZXMgPC0gYXMuZGF0YS5mcmFtZSh0cmFpblssLk4sYnk9Y3VycmVudF9jb2xuYW1lXVtvcmRlcigtTildKVssJ04nXQ0KcGxvdF9seSh4PWFzLmxpc3QoeF92YWx1ZXMpLHk9YXMubGlzdCh5X3ZhbHVlcyksIHR5cGU9ImJhciIpICU+JSBsYXlvdXQodGl0bGU9Y3VycmVudF9jb2xuYW1lKQ0KaTwtaSsxDQpjdXJyZW50X2NvbG5hbWUgPC0gY29sdW1uX25hbWVzX2Zvcl9vcHRpb25fcGxvdHNfc3RyaW5nX2NvZGVzW2ldDQp4X3ZhbHVlcyA8LSBhcy5kYXRhLmZyYW1lKHRyYWluWywuTixieT1jdXJyZW50X2NvbG5hbWVdW29yZGVyKC1OKV0pWyxjdXJyZW50X2NvbG5hbWVdDQp5X3ZhbHVlcyA8LSBhcy5kYXRhLmZyYW1lKHRyYWluWywuTixieT1jdXJyZW50X2NvbG5hbWVdW29yZGVyKC1OKV0pWywnTiddDQpwbG90X2x5KHg9YXMubGlzdCh4X3ZhbHVlcykseT1hcy5saXN0KHlfdmFsdWVzKSwgdHlwZT0iYmFyIikgJT4lIGxheW91dCh0aXRsZT1jdXJyZW50X2NvbG5hbWUpDQppPC1pKzENCmN1cnJlbnRfY29sbmFtZSA8LSBjb2x1bW5fbmFtZXNfZm9yX29wdGlvbl9wbG90c19zdHJpbmdfY29kZXNbaV0NCnhfdmFsdWVzIDwtIGFzLmRhdGEuZnJhbWUodHJhaW5bLC5OLGJ5PWN1cnJlbnRfY29sbmFtZV1bb3JkZXIoLU4pXSlbLGN1cnJlbnRfY29sbmFtZV0NCnlfdmFsdWVzIDwtIGFzLmRhdGEuZnJhbWUodHJhaW5bLC5OLGJ5PWN1cnJlbnRfY29sbmFtZV1bb3JkZXIoLU4pXSlbLCdOJ10NCnBsb3RfbHkoeD1hcy5saXN0KHhfdmFsdWVzKSx5PWFzLmxpc3QoeV92YWx1ZXMpLCB0eXBlPSJiYXIiKSAlPiUgbGF5b3V0KHRpdGxlPWN1cnJlbnRfY29sbmFtZSkNCmk8LWkrMQ0KY3VycmVudF9jb2xuYW1lIDwtIGNvbHVtbl9uYW1lc19mb3Jfb3B0aW9uX3Bsb3RzX3N0cmluZ19jb2Rlc1tpXQ0KeF92YWx1ZXMgPC0gYXMuZGF0YS5mcmFtZSh0cmFpblssLk4sYnk9Y3VycmVudF9jb2xuYW1lXVtvcmRlcigtTildKVssY3VycmVudF9jb2xuYW1lXQ0KeV92YWx1ZXMgPC0gYXMuZGF0YS5mcmFtZSh0cmFpblssLk4sYnk9Y3VycmVudF9jb2xuYW1lXVtvcmRlcigtTildKVssJ04nXQ0KcGxvdF9seSh4PWFzLmxpc3QoeF92YWx1ZXMpLHk9YXMubGlzdCh5X3ZhbHVlcyksIHR5cGU9ImJhciIpICU+JSBsYXlvdXQodGl0bGU9Y3VycmVudF9jb2xuYW1lKQ0KaTwtaSsxDQpjdXJyZW50X2NvbG5hbWUgPC0gY29sdW1uX25hbWVzX2Zvcl9vcHRpb25fcGxvdHNfc3RyaW5nX2NvZGVzW2ldDQp4X3ZhbHVlcyA8LSBhcy5kYXRhLmZyYW1lKHRyYWluWywuTixieT1jdXJyZW50X2NvbG5hbWVdW29yZGVyKC1OKV0pWyxjdXJyZW50X2NvbG5hbWVdDQp5X3ZhbHVlcyA8LSBhcy5kYXRhLmZyYW1lKHRyYWluWywuTixieT1jdXJyZW50X2NvbG5hbWVdW29yZGVyKC1OKV0pWywnTiddDQpwbG90X2x5KHg9YXMubGlzdCh4X3ZhbHVlcykseT1hcy5saXN0KHlfdmFsdWVzKSwgdHlwZT0iYmFyIikgJT4lIGxheW91dCh0aXRsZT1jdXJyZW50X2NvbG5hbWUpDQppPC1pKzENCmN1cnJlbnRfY29sbmFtZSA8LSBjb2x1bW5fbmFtZXNfZm9yX29wdGlvbl9wbG90c19zdHJpbmdfY29kZXNbaV0NCnhfdmFsdWVzIDwtIGFzLmRhdGEuZnJhbWUodHJhaW5bLC5OLGJ5PWN1cnJlbnRfY29sbmFtZV1bb3JkZXIoLU4pXSlbLGN1cnJlbnRfY29sbmFtZV0NCnlfdmFsdWVzIDwtIGFzLmRhdGEuZnJhbWUodHJhaW5bLC5OLGJ5PWN1cnJlbnRfY29sbmFtZV1bb3JkZXIoLU4pXSlbLCdOJ10NCnBsb3RfbHkoeD1hcy5saXN0KHhfdmFsdWVzKSx5PWFzLmxpc3QoeV92YWx1ZXMpLCB0eXBlPSJiYXIiKSAlPiUgbGF5b3V0KHRpdGxlPWN1cnJlbnRfY29sbmFtZSkNCg0KYGBgDQo8YnI+DQoNCiMjIyMgIEZyZXF1ZW5jeSB0YWJsZSBmb3IgdGhlIHJlbWFpbmluZyB2YXJpYWJsZXMgKEFsbCB3aXRoIDAvMSBDb2RpbmcpDQpgYGB7cn0NCmNvbHVtbl9uYW1lc19mb3Jfc3VtbWFyeSA8LSBjb2xuYW1lcyh0cmFpbilbMTE6bGVuZ3RoKGNvbG5hbWVzKHRyYWluKSldDQpWYXJpYWJsZSA8LSBjKCkNCk5fMCA8LSBjKCkNCk5fMSA8LSBjKCkNCk1lYW5fMCA8LSBjKCkNCk1lYW5fMSA8LSBjKCkNCg0KZm9yKGkgaW4gY29sdW1uX25hbWVzX2Zvcl9zdW1tYXJ5KXsNCiAgajwtMQ0KICBWYXJpYWJsZSA8LSBjKFZhcmlhYmxlLGkpDQogIE5fMCA8LSBjKE5fMCx0cmFpblssLk4sYnk9aV1bMV0kTikNCiAgTl8xIDwtIGMoTl8xLHRyYWluWywuTixieT1pXVsyXSROKQ0KICBNZWFuXzAgPC0gYyhNZWFuXzAscm91bmQodHJhaW5bLG1lYW4oeSksYnk9aV1bMV0kVjEsMSkpDQogIE1lYW5fMSA8LSBjKE1lYW5fMSxyb3VuZCh0cmFpblssbWVhbih5KSxieT1pXVsyXSRWMSwxKSkNCiAgaiA8LSBqKzENCn0NCg0KIyBOXzBbaXMubmEoTl8wKT09VFJVRV0gPC0gMA0KIyBOXzFbaXMubmEoTl8xKT09VFJVRV0gPC0gMA0KIyBNZWFuXzBbaXMubmEoTWVhbl8wKT09VFJVRV0gPC0gMA0KIyBNZWFuXzFbaXMubmEoTWVhbl8xKT09VFJVRV0gPC0gMA0KDQpzdW1tYXJ5X3Jlc3VsdHMgPC0gYXMuZGF0YS5mcmFtZShjYmluZChWYXJpYWJsZSxOXzA9Tl8wLE5fMSxNZWFuXzAsTWVhbl8xKSwgc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFKQ0Kc3VtbWFyeV9yZXN1bHRzJE5fMCA8LSBhcy5pbnRlZ2VyKHN1bW1hcnlfcmVzdWx0cyROXzApDQpzdW1tYXJ5X3Jlc3VsdHMkTl8xIDwtIGFzLmludGVnZXIoc3VtbWFyeV9yZXN1bHRzJE5fMSkNCnN1bW1hcnlfcmVzdWx0cyRNZWFuXzAgPC0gYXMubnVtZXJpYyhzdW1tYXJ5X3Jlc3VsdHMkTWVhbl8wKQ0Kc3VtbWFyeV9yZXN1bHRzJE1lYW5fMSA8LSBhcy5udW1lcmljKHN1bW1hcnlfcmVzdWx0cyRNZWFuXzEpDQpzdW1tYXJ5X3Jlc3VsdHMkRGVsdGFfTWVhbiA8LSBzdW1tYXJ5X3Jlc3VsdHMkTWVhbl8xIC0gc3VtbWFyeV9yZXN1bHRzJE1lYW5fMA0Kc3VtbWFyeV9yZXN1bHRzDQpgYGANCjxicj4NCg0KIyMjIyAgQWZ0ZXIgTG9va2luZyBhdCBLYWdnbGUsIGNoZWNrZWQgZm9yIER1cGxpY2F0ZSBmaWxlZHMgLSBBZGRlZCBEYXRhIHRvIEZyZXF1ZW5jeSB0YWJsZSBmb3IgdGhlIHJlbWFpbmluZyB2YXJpYWJsZXMgKEFsbCB3aXRoIDAvMSBDb2RpbmcpLg0KYGBge3J9DQp0cmFpbl8yIDwtIHRyYWluWywgIWR1cGxpY2F0ZWQodCh0cmFpbikpXSAjcmVtb3ZlIGR1cGxpY2F0ZWQgZmllbGRzIC4uLiBmcm9tIHJhZGRhckBLYWdnbGUgPT4gaHR0cHM6Ly93d3cua2FnZ2xlLmNvbS9jL21lcmNlZGVzLWJlbnotZ3JlZW5lci1tYW51ZmFjdHVyaW5nL2Rpc2N1c3Npb24vMzQwMDYNCmR1cGxpY2F0ZV9jb2x1bW4gPC0gYygpDQpmb3IoaSBpbiBkdXBsaWNhdGVkKHQodHJhaW4pKSl7DQogIGR1cGxpY2F0ZV9jb2x1bW4gPC0gYyhkdXBsaWNhdGVfY29sdW1uLGkpDQp9DQoNCmR1cGxpY2F0ZV9jb2x1bW4gPC0gZHVwbGljYXRlX2NvbHVtblsxMTpsZW5ndGgoZHVwbGljYXRlX2NvbHVtbildDQoNCnN1bW1hcnlfcmVzdWx0cyRkdXBsaWNhdGVfY29sdW1uIDwtIGR1cGxpY2F0ZV9jb2x1bW4NCg0Kc3VtbWFyeV9yZXN1bHRzX2RlY3JlYXNpbmcgPC0gc3VtbWFyeV9yZXN1bHRzW29yZGVyKC1OXzEpLF0NCg0Kc3VtbWFyeV9yZXN1bHRzDQpzdW1tYXJ5X3Jlc3VsdHNfZGVjcmVhc2luZw0KYGBgDQo8YnI+DQoNCiMjIyMgIExpc3Qgb2YgRHVwbGljYXRlZCBDb2x1bW5zDQpgYGB7cn0NCnN1bW1hcnlfcmVzdWx0c19kZWNyZWFzaW5nWyhzdW1tYXJ5X3Jlc3VsdHNfZGVjcmVhc2luZyRkdXBsaWNhdGVfY29sdW1uPT1UUlVFKSxdDQpgYGANCjxicj4NCg0KIyMjIyAgQ3JlYXRlIFJlc3BvbnNlIFBsb3RzIGZvciBhbGwgQmluYXJ5IE9wdGlvbnMNCmBgYHtyfQ0KdHJhaW5fZGYgPC0gYXMuZGF0YS5mcmFtZSh0cmFpbikNCnRyYWluX2RmX2RlZHVwZSA8LSB0cmFpbl9kZlssICFkdXBsaWNhdGVkKHQodHJhaW5fZGYpKV0gI3JlbW92ZSBkdXBsaWNhdGVkIGZpZWxkcyAuLi4gZnJvbSByYWRkYXJAS2FnZ2xlID0+IGh0dHBzOi8vd3d3LmthZ2dsZS5jb20vYy9tZXJjZWRlcy1iZW56LWdyZWVuZXItbWFudWZhY3R1cmluZy9kaXNjdXNzaW9uLzM0MDA2DQoNCmNvbHVtbl9uYW1lc19mb3JfYmluYXJ5X3Bsb3RzIDwtIGNvbG5hbWVzKHRyYWluX2RmX2RlZHVwZSlbMTE6bGVuZ3RoKGNvbG5hbWVzKHRyYWluX2RmX2RlZHVwZSkpXQ0KDQpsIDwtIGh0bWx0b29sczo6dGFnTGlzdCgpDQpmb3IgKGkgaW4gMTpsZW5ndGgoY29sdW1uX25hbWVzX2Zvcl9iaW5hcnlfcGxvdHMpKSB7ICAjDQogIHBsb3RfdmFyaWFibGUgPC0gY29sdW1uX25hbWVzX2Zvcl9iaW5hcnlfcGxvdHNbaV0NCiAgbFtbaV1dIDwtIHBsb3RfbHkoeD10cmFpbl9kZl9kZWR1cGVbLHBsb3RfdmFyaWFibGVdLHk9dHJhaW5fZGZfZGVkdXBlJHksIHR5cGU9ImJveCIsIGhlaWdodD0zMDApICU+JSBsYXlvdXQodGl0bGU9cGxvdF92YXJpYWJsZSkNCn0NCmwNCg0KYGBgDQo8YnI+DQoNCg==